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Abstract 

A network of agents cooperate on a given area. Time evolution of 
their power is described within a set of nonlinear equations. The limi- 
tation of resources is introduced via the Verhulst term, equivalent to a 
global coupling. Each agent is fed by some other agents from his neigh- 
borhood. Two subsequent stages of the time evolution can be observed. 
Initially, the richness of everybody increases distinctly, but its distribu- 
tion becomes wide. After some transient time, however, resources are 
exhausted. Richness of some agents falls to zero and they are eliminated. 
Cooperation becomes less effective, what leads to subsequent falls. Fi- 
nally, small percent of agents survive in a steady state. We investigate, 
how the cooperation influences the rate of surviving. 

1 Introduction 

Physicists are happy. Their resources — amount of problems to solve — is 
infinite. It is not so, however, in almost all other professions; the numbers of 
car buyers, voters, butterflies to catch and girls to kiss for the first time are 
limited. If one takes everything, others have to look for another hobby or job. 
On the other hand, cooperation is a fingerprint of modern society. All what we 
get except from our domestic gardens, we get from other people. The aim of this 
work is to present numerical results of a set of equations, where the above facts 
are built in as assumptions. The only way to increase one's power or richness 
or speed of getting resources, which are treated here as synonymous, is to profit 
work of somebody else. Then, each agent i has M neighbours who feed 

its with given speeds Jy. The limitation of resources is taken into account as 
the global coupling via a nonlinear Verhulst-like term. These are first two r.h.s. 
terms of our basic equations (Q), and these terms do not depend on i. Third 
term is to include some action of an i-th agent. Namely, it selects from his 
neighbours the one most sensitive for this action, and enhances feeding from 
this particular neighbour. 

The problem belongs to the large class of models designed to describe a 
competition for resources Q. Most frequently, however, some kind of dynamic 
equilibrium is considered, and only a few authors are interested in an ultimate 
catastrophe. Example giving, although statistical physics provides tools for 
analysing stock market |2], the term "bankruptcy" is absent in physical journals 
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(see [|| and for exceptions). It seems worthwhile to take a glance on a collision 
of expanding society with the boundary of limited resources. 

2 Model 

Time evolution of the individual agent power Pj is described by differential 
equation: 



' WAh -A 1 -^J i ,P J -(i)r J -(t)-A 2 



dt 

3=1 



1 N 



(1) 



+A 3 ■ max [di,-Pj(t)r,-(i)] 



Coupling constant J^- (speed of feeding) and d^ (sensitivity) are random posi- 
tive reals normalized to unity and fixed during simulation, while Ai, A 2 and A3 
describe intensities of the three terms. The constants r.i are equal to one for 
active agents and to zero in other case. At the beginning, all agents are active, 
but once Pi is negative for given i, ri is switched to zero and the ?-th agent is 
eliminated from the game. We deal with a set of differential equations which are 
piecewisely continuous. At the moments of time t* when Pi(t*) — for any up- 
to-now-active agent i, the equations are switched from one analytical solution to 
another one. In this sense, the formalism is equivalent to a coupled map lattice, 
but the number of equations changes in time. Similar approach was applied 
already in j|. Note however that in our case, the maps are to be integrated nu- 
merically. We consider two different neighbourhoods: (i) each individual has M 
randomly chosen neighbours, or (ii) each individual has M nearest geometrical 
neighbours. In the latter case, the agents form a one-dimensional chain with 
periodic boundary conditions. Note that the formed network of agents is a con- 
tinuous analogue to the Kauffman model Q , designed for a simulation of genetic 
systems. We solve (Q) numerically and check how the kind of neighbourhood, 
the number of directly interacting agents M and the strengths of interactions 
influence (i) time evolution of the average power (ii) and number of active agents 
after very long times. The simulation is started with N = 10 4 active agents 
(n = 1 for i = 1, • • • , N ), each with randomly chosen initial power Pi{t = 0). 
The simulation takes Ai ter = 10 4 time steps, each 10 -3 long, what guarantees 
the numerical stability. 



3 Results and discussion 

Two subsequent stages of the time evolution can be observed independently of 
kind of neighbourhood (Fig. |l|). Initially, the average power increases distinctly, 
but its distribution becomes wide. After some transient time, however, resources 
are exhausted. Richness of some agents falls to zero and they are eliminated. 
Cooperation becomes less effective, what leads to subsequent falls. Finally, 
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Figure 1: Time evolution of (a) the average agent power (P) expressed in arbi- 
trary units and (b) subsequent percentage n of active agents for random (doted 
line) and geometrical (solid line) neighbours. 
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Table 1: The percentage of the society's successes for geometrical neighbourhood 
and various sets of (Ai, A3). M = 4, A 2 = 1, No = 10 4 , Mtcr = 10 4 , iV run = 1000. 



small percent of agents survive in a steady state. Distinct differences are found 
between the cases of geometrical and random neighbourhood (Fig. |l|). For 
the case of random neighbours, usually only one agent survives, i — ult. This 
simplifies Eq. (|l|) to 

^§^ = -A 2 -P ult W (2) 

what gives subsequent power decrease of the last agent as 1/t. If two or more 
cooperating agents survive, we are faced with a set of nonlinear equations. An- 
alytically, a simplified case can be considered, when Jy = J, dy = d. Then, the 
average power tends to a positive stable fixed point. Numerical results suggest, 
that this is the rule also in the general case (Fig. [j]). The marked difference 
between the results for random and geometrical neighborhoods is an illustration 
of the old truth do ut des. In other words, it is better to help friends which 
can reward than to people randomly selected in the street. A group of agents 
feeding each other can survive, if they spend resources moderately enough. A 
kind of an equilibrium with a given environment seems to be possible for the 
geometrical neighborhood, as long as the Verhulst term is compensated by the 
remaining terms in Eq. ([j]). 

The results given above refer to the case of the geometrical neighborhood. 
Let us define as a success of a society the case when more than one cooperating 
agents survive. Then, Tab. [I] gives the percentage of successes for various sets 
of (Ai, A 3 ). It is astonishing (at least for us) that the roles of these two terms 
(first and third) in the r.h.s. of Eq. |l]) differ so much. Actually, both these 
terms are designed to increase the power of a given agent. The difference is 
only that the third term is a kind of local optimization, while the first one 
is automatic. It seems that this kind of dynamic reaction of a given agent is 
particularly relevant at the border of "death" (or "bankruptcy" or so), when 
Pi is close to zero. Average power in asymptotic time, presented in Tab. 
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Table 2: Average power (P) (in arbitrary units) after Ai tcr = 10 4 time steps 
averaged over N run = 100 independent runs for geometrical neighbourhood and 
various (Ai,A 3 ). M = 4, A 2 = 1, N = 10 4 . 



does not show this effect. There, both terms act in the same way and can be 
mutually replaced to get approximately the same result. 

4 Conclusions 

We are somewhat surprised with the fury of crisis, which can be observed in 
Fig. [I]. Simultaneously, distinct but continuous progress of the average power is 
substituted by wild oscillations, and the number of agents is strongly reduced. 
The effect arises abruptly, without any preceding warnings in the curve shape. 
Soon, unavoidable elimination of almost all agents is observed for all applied 
sets of input parameters. We imagine the model as a parabolic description of 
the process of breaking of bonds and a destruction of a complex system. 

We feel that in this text, full of analogies, we are at the border of abuse of 
language. We apologize for that. However, we have found that it is particularly 
difficult to describe a transient process in precise terms of statistical physics, 
most of them designed for stationary processes. Our process could be treated as 
stationary, if we allow the population of agents to be slowly reproduced, maybe 
with retaining some fruitful information, and the whole story is repeated many 
times. This kind of investigation could link to a Penna-like model M and to 
some kind of self-organization Q of the structure and strength of bonds between 
agents. 
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